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Abstract: We demonstrate that due to strong modal interactions through 
cross-gain saturation, the onset of a new lasing mode can switch off an exist¬ 
ing mode via a negative power slope. In this process of interaction-induced 
mode switching (IMS) the two involved modes maintain their identities, 
i.e. they do not change their spatial field patterns or lasing frequencies. 
For a fixed pump profile, a simple analytic criterion for the occurrence 
of IMS is given in terms of their self- and cross-interaction coefficients 
and non-interacting thresholds, which is verified for the example of a 
two-dimensional microdisk laser. When the spatial pump profile is varied as 
the pump power is increased, IMS can be induced even when it would not 
occur with a fixed pump profile, as we show for two coupled laser cavities. 
Our findings apply to steady-state lasing and are hence different from 
dynamical mode switching or hopping. IMS may have potential applications 
in robust and flexible all-optical switching. 
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1. Introduction 


Our understanding of lasers is largely based on the semiclassical laser theory, which describes 
successfully modal thresholds, output power, spatial mode patterns inside and outside of the 
laser cavity, as well as dynamical effects such as relaxation oscillations and mode, phase, and 
frequency locking W- With the advent of microlasers their non-Hermitian nature and 
strong modal interactions not only introduced new phenomena (see, for example, |j^|^) but 
also imposed new challenges for the understanding of these novel lasers. A new framework 
based on the semiclassical laser theory, known as the Steady-state Ab-initio Laser Theory 
(SALT) |[6|{10|, addresses these challenges satisfactorily when the level inversion in the gain 
medium can be treated as stationary. SALT has been shown to give excellent agreement with 
much more computationally intensive finite-difference-time-domain (FDTD) simulations of the 
lasing equations |TT]-T^, and it also accurately predicts recent experimental observations of 
exceptional points in lasers @0- These advantages make SALT an ideal tool to study the 
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Fig. 1. Schematics showing modal intensities as a function of the pump strength for (a) 
interaction-induced mode switching and (c) typical two-mode lasing. (h,d) The correspond¬ 
ing lasing spectra at pump strengths marked hy I and II in (a,b). 


effects of non-linear modal interactions in the steady-state, in regimes not previously studied 
and in some cases not accessible with other approaches. 

In this report, we reveal a surprising phenomenon due to such modal interactions: the onset of 
one lasing mode can switch off another mode which is already lasing [Figs.[T|a) and[T|b)]; we 
refer to this phenomenon as interaction-induced mode switching (IMS). IMS is very different 
from typical multimode lasing, where the modal intensity of each mode increases linearly with 
the pump strength and the power slopes of lasing modes are only slightly reduced at the onset of 
a new mode [Figs.[TJc) and[TJd)]. IMS is a robust phenomenon and can in principle occur in 
any type of laser, including standard systems such as a microdisk and coupled one-dimensional 
cavities. We note that IMS is a steady-state effect and different from dynamical mode switching 
and related phenomena | T5p^ . It is also independent of laser bistability 117 ■ which can also 
lead to various mode switching scenarios. In addition, IMS differs qualitatively from seemingly 
similar effects found in our previous work I00, as we will explain in detail below. 

In the next section we hrst review the basic ansatz and results of SALT and its “single-pole” 
approximation (SPA-SALT), where each lasing mode can be well approximated by a linear 
mode, the constant-flux (CF) state QSIID- We then present an example of IMS in a uniformly 
pumped two-dimensional (2D) microdisk laser. The presence of IMS in this example can be 
predicted and understood via a simple analytic criterion obtained from SPA-SALT, involving 
only a few parameters, i.e., the self- and cross-gain saturation coefficients of the two modes 
(“interaction coefficients”) and the thresholds of the modes in the absence of gain saturation 
(“non-interacting thresholds”). The prediction of IMS from SPA-SALT is conhrmed by the use 
of the full SALT computation. Next, we show that IMS can also be facilitated by having a spatial 
pump prohle which evolves as the overall pump power is increased. An example is given using 
two coupled one-dimensional (ID) slab cavities, in which IMS does not occur when the cavities 
are pumped uniformly. With different choices of evolving the spatial pump prohle, we show that 


























either of two modes can undergo an IMS, i.e. mode 1 be switched off by mode 2 or vice-versa, 
indicating the possibility of flexible control of the mode spectrum by using IMS and variable 
pump profile. A criterion for IMS in this case is again given using SPA-SALT, which shows an 
invariant property of IMS and its less restrictive requirement on the interaction coefficients. 


2. SALT and SPA-SALT 


SALT makes a steady-state ansatz and reduces the time-dependent coupled lasing equations 
for the electric field, level populations and polarization of the gain medium, to a set of non¬ 
linear time-independent wave equations for the lasing modes and their frequencies |[6 ■ 121. In 
the case of multimode lasing it takes advantage of a separation of time scales between the 
slow population dynamics and the beat frequencies of the modes to predict a stable multimode 
stationary state, an approach which is typically valid only for microlasers |T^^|23[. While 


SALT neglects the time-dependent effects of the four-wave mixing non-linearity, it treats the 
time-independent effects of gain saturation and spatial hole-burning to infinite order and hence is 
found to be accurate far above threshold in steady-state lasing. In the discussion below, diffusion 
of excitations in the gain medium is assumed absent; recent work has shown how to include this 
effect in the SALT framework pO) . 

SALT differs from the standard modal description of semiclassical laser theory in a few 
important aspects: SALT describes the openness of the lasing medium rigorously in terms of 
an outgoing boundary condition, and, as noted, the static modal interaction is treated to infinite 
order and not truncated at the cubic level. The first method developed for solving the SALT 
equations and the one most useful in the current work is through the expansion of a lasing 
mode 'P^(r;Z)o) at a given pump strength Dq in an optimal basis of CF states {up{r, COp)}, i.e. 

{r,Do) = Lp^ {Do)up{r,C0p). Each CF state is purely outgoing and parametrized by the 
lasing frequency cOp, and the full set is biorthogonal and complete; hence an arbitrary lasing mode 


can be represented in this fashion. We define the modal intensity by Ip {Do) = 


\Do)\^ 


and its power slope, Sp = dip {Do)/dDo, is the key quantity of interest in the study of IMS. 

In the previous expressions, the lasing frequency (Op is a real quantity to be solved together 
with the mode profile '¥p{r). The pump strength is represented as a dimensionless quantity Do 
proportional to the level inversion in the gain medium, and its spatial profile is assumed given 
by f{r) and normalized by /cavity“ /cavitybiorthogonality relation of the CF 
states can be shown to be 0 


/ dru*p{r)uci{r) =V5pci 
J cavity 


( 1 ) 


when the refractive index n is uniform inside the cavity. Up{r) is the biorthogonal partner of Up{r) 
and “ * ” denotes complex conjugation (for a lossless cavity without degeneracy, the partner is 
simply Up{r) itself). In the next section we consider degenerate CF states, the traveling-wave 
(clockwise (CW) or counterclockwise (CCW)) modes in a microdisk, i.e.. 
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In this case Up{r) is then given by 
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Here 0 is the azimuthal angle, is the Bessel function of azimuthal quantum number m, D.p is 
the complex CF frequency, c is the speed of light in vacuum, and R is the radius of the microdisk 
cavity. 
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When the Q-factors of the lasing modes are relatively high, each lasing mode can be well 
approximated by a single CF state, i.e., there is only one significant in the expansion of 
'Fu, even when well above the laser threshold and when the pump profile f{r) is non-uniform 
. This property enables us to invoke SPA-SALT, i.e., 'P^(r;Do): 
tfiiDo) « This approximation neglects the power-dependence of the lasing frequencies 

and the change in the spatial mode profiles above threshold, but has the great advantage of 
reducing the non-linear coupled SALT equations to a set of (constrained) linear equations for 
the modal intensities Q. For example, in the two-mode regime, the modal intensity Ii and I 2 are 
determined by 
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(4) 


with the constraint that /i ,/2 > 0. Here F^ = -f — (Oaf'] < 1 is the Lorentzian gain 

factor for mode /t, where (Oa is the atomic transition frequency and Y± is the longitudinal 
relaxation rate of the gain medium. The self-interaction coefficients and the cross¬ 

interaction coefficients Xi2,X2i are given by 


Xpij — 
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We note that the product nClp is almost real for a relatively high-Q mode, even when the passive 
cavity refractive index n has a small imaginary part due to the inclusion of material absorption. 
Therefore, using the expressions 0 and ([^, the cross-interaction coefficients are approximately 
equal, i. e., 

'nQ,ir\ / nQ,2r'' 
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J cavity 

in the microdisk laser to be discussed in the next section. 

Since we are using the microdisk laser as an example, we need to discuss the conditions 
for the validity of SALT (and SPA-SALT) in more detail. The rates which determine the laser 
population dynamics (in the simplest case of the Maxwell-Bloch model of two-level gain atoms) 
are described by a single parameter, denoted as Y\\■ Multimode SALT is approximately valid when 
7 x, A ^ Y\\ where A is the typical beat frequency between lasing modes (of order 

the free spectralrange). For the degenerate CW and CCW modes of a microdisk this condition is 
not satisfied, and we will assume that, as is typical in ring lasers, one sense of rotation is chosen 
by small perturbations, and SALT can then be used for the resulting non-degenerate lasing 
modes. A more detailed treatment of the degenerate case can be found, for example, in | [2^ , but 
it is not relevant for the current discussion. We also note that the electric field and level inversion 
are measured in their respective namral units /'^8 dc = hY_i_/4ng^, where g is 

the dipole matrix element between the lasing transition levels. 

The key quantity of interest in the discussion of IMS in the two-mode regime is the power 
slope of the first mode, as we have shown schematically in Fig.[TJa). Within SPA-SALT all lasing 
modes have a linear variation with pump power if the spatial pump profile is fixed, featuring 
a slope which only changes when new modes turn on or off. For the two mode case, it is 
straightforward to show that the power slope of the first mode is Si = l/[ri;KiiDQ'^] after it 
turns on at Dq = in which Xii plays the role of mode volume; this slope changes to 
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once the second mode starts lasing at 
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The subscript “int” here indicates that Dg is the threshold of mode 2 in the presence of modal 

interaction with mode 1. In contrast, in Eqs. Q and ([^ is the threshold of mode 2 when the 
modal interaction is neglected, which is determined by the openness of the laser cavity and the 
material absorptions within, as well as the Lorentzian gain factor We will refer to 

as the non-interacting threshold of mode . Since there is only one lasing mode at the first 


threshold, Og\^^^ = Og^ by dehnition. 

Within SPA-SALT the full lasing behavior is captured by the set of interacting thresholds 
and power slopes between thresholds; these are determined solely by the non-interacting lasing 
frequencies and thresholds, the gain curve parameters and the interactions coefficients. By 
definition a mode which has just turned on has a positive power slope; IMS occurs when at some 
interacting threshold the power slope of a previously lasing mode becomes negative, e.g., < 0 

in the two-mode case above. This criterion is completely general and can be applied to any laser 
system for which SPA-SALT is a good approximation. 

The SPA-SALT equations Q and their multimode form may look similar to the Haken- 
Sauermann equations |22||24| and their generalizations (e.g., |25j), but the latter employ the 
cubic approximation for the modal interactions and lead to unphysical mode behaviors for pump 
strength ^ 20% above the first threshold Q. A detailed discussion of these points can be found 
in @[ID, where SALT is also compared to the work by Mandel and coworkers p^[27[. 


3. IMS with a fixed pump profile 

In this section we give an example of IMS in a 2D microdisk laser. We consider uniform pumping 
(f(r) = 1) and transverse magnetic modes, with 'P^(r) representing the out-of-plane electric 
field in the steady state. We restrict our treatment to CCW traveling waves solutions for the 
reasons discussed previously. 

Using the full SALT, we find that a whispering gallery mode of azimuthal quantum number 
mi =8 starts lasing when the pump strength is higher than its threshold Dg^^ = 0.075 (see Fig.|^. 
Its modal intensity increases linearly with the pump strength, until a second mode of m 2 = 7 
turns on at Dg = 0.17, whose non-interacting threshold is at Dq = 0.076. Beyond Dg an 
anomalous behavior occurs: unlike in a typical laser where the onset of a new mode just slightly 
slows down the increase of the intensity of the the existing mode(s) [Fig.[2c)], the power slope 
of mode 1 actually becomes negative, and it stops lasing at = 0.19. Note that we have used 
here full SALT, for which the intensity variation with pump need not be linear Q, but in fact it 
is quite linear as predicted by the SPA-SALT approximation. This phenomenon is confirmed 
with good agreement by directly solving the differential form of SALT numerically | [T3) , without 
the expansion in the CF basis (not shown). 

To understand the origin of this mode switching behavior we analyze in more detail the 
criterion for IMS in the two mode case derived from the SPA-SALT approximation, that 
§1 < 0. Let us first discuss the typical behavior in a two-mode laser shown schematically 
in Fig. [TJc). One usually finds the cross-interaction coefficients to be much smaller than the self¬ 
interaction coefficients, i.e., Xi2,X2i ^ Zi 1 jZ 22 - Therefore, for two lasing modes with similar 
non-interacting thresholds, we immediately see that from Fq. (j^, which indicates that 

the onset of the second lasing mode slightly reduces the power slope of the existing mode, and 
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Fig. 2. Interaction-induced mode switching in a microdisk laser. Lower panel: Black solid 
and red dashed lines show the modal intensities of the first two modes. They are whispering- 
gallery modes with azimuthal quantum numbers mj = 8, m 2 = 7 and lasing frequencies 
a\R/c = 5.37, OiiRlc = 4.81, respectively. Upper panel: False color plots of the real part 
of the electric field in these two modes, where red and blue indicate positive and negative 
values. Their radial profiles are shown as the inset in the lower panel. The parameters used 
are: atomic transition frequency cOaR/c = 4.83, longitudinal relaxation rate y^R/c = 1 and 
refractive index n = 2 + 0.01 1 . 


§1 remains positive. Consequently, the IMS behavior shown in Fig.|^must be due to a strong 
cross-interaction. 

We now derive the relevant criterion for IMS in terms of the self- and cross-interaction 
coefficients. We first note that by definition ^ ^ 0 ^^’ i-®-’ second mode turns on at 

a higher pump strength than the first mode, with or without modal interactions. Therefore, we 
find 




Xu D, 
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from Eq. (|^. With this observation, we derive from Eq. Q the following criterion for a negative 
5i: 
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When satisfied, mode 1 is switched off at 
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This expression is consistent with Eq. (|8 1 , which leads to > T>o7nt’ 1® ’ termination of 
mode 1 happens after the onset of mocM 2. 

The criterion ( [T0| ) requires a cross-interaction coefficient larger than the self-interaction 
coefficient of the “killer” mode 2 and smaller than that of the “victim” mode 1: 


,( 2 ) 


X22<X2i^Xi2<Xn- (12) 

To find the circumstances in which these inequalities can hold, we note that in the definition 
(|^ of the interaction coefficient Xpq, the intensity profile of mode p (q) can be treated as the 
weighting function for mode q (p). To have a large Xpq, the best weighting function would be a 
delta function at the peak(s) of mode q {p). Therefore, for the criterion ( p^ to hold, the intensity 
profile of mode 1 needs to be peaked very close to the peak(s) of mode 2, and it should be more 
localized. Eor the two whispering-gallery modes shown in Eig.|^ the first mode {m\ — 8) indeed 











Fig. 3. (a,c) Same as the modal intensity plot shown in Fig.j^hut with COaR/c = 0.50 (a) and 
(OaRjc = 0.48 (c). Single-mode lasing is observed and mode switching does not occur in 
either case, even when the pump power is high above threshold. (b,d) The modal gain of 
these two modes in (a,c), confirming that the m 2 = 7 mode in (a) and the mi = 8 mode in (c) 
do not lase. Modal gain is clamped at 1 for lasing modes and stays below 1 for non-lasing 
modes Q. 


has a slightly narrower peak in the radial direction that overlaps with the peak of the second mode 
(see the inset). As a result, we find X 21 /X 11 — 0.977, X 22 /Z 12 = 0.979, and Eq. (12 1 is satisfied. 
In addition, we find that the complete criterion (10 1 also holds, because /D^’= 0.981. 
Due to the small differences between the ratios X 21 /Z 11 , XiilXii, and D\ 


(1) /n(2) 


ID\ 


which 


appear in the denominators of Eqs. ® and ( [TT] ), the quantitative values of the onset threshold 

'—' (2) 

of mode 2 and the termination threshold of mode 1 as given by SPA-SALT (DJj ’ 


0.41, 
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( 2 ) . 
l.off ' 


: 0.70) are susceptible to the inaccuracy introduced by the single-pole approximation and 


("T] z. 

do not agree very well with their acmal values {D^ « 0.17, « 0.19). Nevertheless, the 

predictive power of Eq. ( [T0| can be confirmed by showing the absence of IMS when Eq. ( [T0| i is 
violated by some margin. 

We can most conveniently violate the SPA-SALT criterion for IMS by varying the last 
ratio of Eq. ( |10| l in our numerical calculation; we achieve this by moving the gain center (Oa- 
This approach changes the non-interacting thresholds of these two modes but barely affects 
their lasing frequencies and mode patterns. Hence j’s and the first two ratios in Eq. (10 1 can 
be treated as unchanged. We first move the gain center closer to the mi = 8 mode, which 
reduces its non-interacting threshold while increasing that of the m 2 = 7 mode. Eventually 
the last ratio in Eq. ( |10| l becomes smaller than the middle one and Eq. ( |10| l breaks down. 
In Fig. ^a) we show such a case with (Oa = 5.0, where = 0.0663,Dg = 0.0792, and 

= 0.837 < Xii!Xw^Xii!X\2- Indeed we find that IMS does not take place, and the 
m 2 = 7 mode is suppressed even when the pump strength is 10 times its non-interacting threshold. 
Next we move the gain center closer to the m 2 = 7 mode, which becomes the first one to lase. 
Thus the necessary condition •EH breaks down, since now the mode index 1 and 2 are exchanged. 
Therefore, the criterion ( [T0) i also breaks down, and again we find that IMS does not take place as 
a result. In Eig.j^c) we show such a case with (i>a = 4.80, where the non-interacting thresholds 
for the m = l modes are 0.0762, 0.0767, respectively. As a conclusion, we find that Eq. 
from the SPA-SALT derivation gives a reliable criterion for the occurrence of IMS. 














We also note that in Fig. we have chosen the gain center to be in close vicinity of the 
frequency of the second mode. This choice was made such that the non-interacting threshold 
of the second mode, which has a lower Q-factor, is not too high to violate the last inequality of 
Eq. ([T§. However, IMS in general does not require the gain center to be closer to the second 
mode than the hrst mode. 

Finally, we address in more detail the difference between the IMS phenomenon studied here 
and the “mode-killing” effect found previously in 2D quasi-ballistic random lasers (QBRL) Q. 
Because QBRL have very low Q, the lasing modes change considerably in their spatial mode 
patterns above threshold and their lasing frequencies have strong power dependence, both of 
which are not included in SPA-SALT. It was found that when two lasing frequencies began to 
cross as a function of power, one of the solutions would reduce in intensity to zero, but the 
remaining mode would change dramatically its spatial pattern and would evolve into a fairly 
equal mixture of the two modes before the frequency crossing. This “mode-mixing” scenario 
requires very low Q and does not occur even in more disordered diffusive random lasers. It 
cannot play a role in the current example for two reasons. First, with uniform pumping the 
microdisk laser conserves the azimuthal quantum number m, so our two interacting modes of 
different ins do not mix spatially through the non-linear interaction; they only affect each other 
through gain saturation. Second, the frequencies of these two modes are well separated and 
do not change with increasing pump. In fact we would not expect IMS to occur in general in 
random lasers, because we require quite large cross-interaction coefficients, which are not likely 
to occur between two modes of random lasers. 

4. IMS with an evolving pump profile 

In this section we show that when IMS does not occur for a fixed pump prohle, it can be 
facilitated by a pump prohle that evolves with the pump strength; IMS can then be induced as a 
result of the combined actions of a linear effect (i.e., change of non-interacting thresholds with 
the varied pump prohle) and the non-linear modal interaction. 

Conceptually IMS is less surprising in this case. First one chooses a system in which different 
modes are to some extent spatially distinct, so that by changing the spatial prohle of the 
pump one can favor one mode over another with similar Q value. Then one simply dehnes a 
“pump trajectory” Q which initially favors one mode and then evolves so as to favor the other. 
Eventually the effective gain of the second mode becomes so high that it forces the hrst mode 
below threshold by gain saturation m- 

However, this intuitive picture does not always hold, because there is no guarantee that the 
effective gain of the hrst mode will reduce eventually as the overall pump power increases 
(see the quantitative discussion below Eq. ( [T6) l near the end of this section). More importantly, 
evolving the pump prohle is an important generalization because one relies much less on hne 
tuning of parameters of the cavity and gain medium to generate IMS, which then occurs more 
robustly. This is exemplihed in Fig.|^using two coupled ID laser cavities with slightly different 
lengths (Li < Li) and independent pump controls (Di ,^ 2 ). The gap between the two cavities 
(Vk) is much shorter than their lengths, and there exist pairs of lasing modes that have similar 
frequencies but reside asymmetrically in these two cavities. 

One such pair (mode 1 and 2) are shown in Eig.|^a). Mode 1 has a stronger intensity in the 
left cavity, while mode 2 does the opposite. If the two coupled cavities are pumped uniformly, we 
hnd that mode 1 turns on hrst at Do = 0.313, and its intensity continues to increase after mode 2 
turns on at Do = 0.349 [see Eig.|^c); inset]. The absence of IMS in this case is well captured by 
the criterion ( [l^ , which does not hold since X 22 /Xi 2 = 1-64 is larger than = 0.953. 

To generate IMS for this system and these two modes, we evolve the pump prohle as we 
increase its average strength, by uniformly pumping each cavity but with different weights. 
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Fig. 4. Interaction-induced mode switching facilitated by an evolving pumping in two 
coupled ID cavities of index n = 3. Their lengths are Li = 4.2/ira, L2 = S.Oflm and the 
gap width between them is IV = O.&flm. L = L\ + Li + W = Wflm. (a) Intensity profiles 
of two modes with wavelength X\ = 4.15/lm, A 2 = 4.28/lm. Shaded areas show the two 
cavities, (b) Their non-interacting thresholds (solid and dashed lines) versus the pump ratio 
D 2 /D 1 . The average pump strength Dq along the pump trajectory used in (c) is shown as 
the dash-dotted line, and the arrow indicates the direction of increasing D\. The triangle 
and diamond on it mark the onset of mode 2 and the termination of mode 1 . The gain curve 
is centered at (HaLjc = 15 = 4.19/rm) and its width is y^L/c = 1. (c) IMS of mode 1 

when the pump strength in the right cavity is increased while the pump strength in the left 
cavity is fixed at Dj = 0.318. Inset: Mode switching does not occur for uniform pumping 
(D\ = D 2 ). (d) IMS of mode 2 when the pump strength in the left cavity is increased while 
the pump strength in the right cavity is fixed at D 2 = 0.344. Good agreement between SALT 
(lines) and FDTD simulations (triangles) are shown. y\\/y± = 5 x 10^^ is used in FDTD 
simulations. 


Specifically, the pump profile f{x) is normalized by f{x)dx = f\Li +/ 2 i 2 = +L 2 , and 

the average pump strength is given by Dq = {D\Li +D 2 L 2 ) / +L 2 ), with £> 1 ( 2 ) = £^ 0 / 1 ( 2 ) by 
definition. Since mode 1 is stronger in the left cavity, its threshold increases if we pump the 
right cavity preferably [see Fig.Qb)]; the opposite holds for mode 2, because it has a stronger 
intensity in the right cavity. For a fixed ratio / 2//1 < 1.11 (including the uniform pumping case), 
mode 1 is the first mode to lase when the pump strength is increased. 

Starting from uniform pumping and Di — D 2 = 0.318 (grey dot in the inset of Fig.j^c)), we 
increase D 2 while keeping Di fixed. The modal intensity of mode 1 continues to increase until 
the onset of mode 2 at D 2 = 0.342 (/ 2//1 = 107), after which it reduces with D 2 [Fig.|^c)]. 
Mode 1 is switched off eventually at D 2 = 0.385 (/ 2//1 = 1.21). It is clear that the non-linear 
modal interaction still plays a crucial role in this mode switching scenario; mode 1 would have 
continued to lase above D 2 = 0.385 if not interacting with mode 2, because the average pump 
strength is still above its non-interacting threshold [Fig.j^b)]. This behavior is distinct from the 
self-termination of a lasing mode due to an exceptional point as reported in ||5 14 the latter 
is a linear effect and does not depend on the non-linear modal interactions. 

Interestingly, with this additional control knob, the IMS scenario can be interchanged if 
a different initial pump configuration is chosen. In Fig. |^d) we start with D\ — 0.287 and 
D 2 — 0.344 (D 2 /£)i = 1.2). Only mode 2 lases since the pump strength is just above its non¬ 
interacting threshold but below that of mode 1. As we increase Di while keeping D 2 fixed, the 
intensity of mode 2 increases linearly until D\ > 0.306 (D 2 /£)i < 1.12), beyond which mode 1 















turns on and switches off mode 2 at Di = 0.331 {Di/Di = 1.04). 

For this one-dimensional example it is relatively easy to provide a test of our results indepen¬ 
dent of SALT by using time integration of the full Maxwell-Bloch equations iTTJT^. As shown 
in Fig.|^d), such calculation also predicts IMS and agrees with SALT, finding less than a 1% 
difference between the values of the pump strength at which mode 1 is killed. This difference 
is within the typical range (a few percent) when comparing SALT and FDTD results 111 T?). 
We also note that our coupled-cavity laser does not show signs of bi-stability |[T7}|^, i.e., the 
frequencies and intensities of the lasing modes are unique for a given pump combination Di and 
D 2 , even though it can be reached by evolving different initial pump configurations. One such 
comparison is shown by the grey dots in Fig.j^c) and its inset, which give the same intensity of 
mode 1 at Di = D 2 = 0.318. Another comparison is shown by the black dots in Fig.Qd) and the 
inset of Fig.j^c), which give the same intensity of mode 1 at Di = D 2 = 0.344. The absence of 
bistability here is due to the conditions 7 x, A ^ we consider, with which a stationary inversion 
can be assumed in the two-mode regime |j7 11 13 it is exact in the single mode regime. 

Besides the flexibility of exchanging the killer mode and the victim mode, IMS is also robust 
and exists in a wide parameter space when evolving the pump profile. For example, IMS shown 
in Fig.|^c) does not rely on the specific value of the fixed pump strength Di; it always occurs as 
D 2 is increased, provided that the mode 1 still lases first. Surprisingly, we find that the modal 
intensity Ii along these different pump trajectories are all invariant in shape as a function of D 2 , 
as shown in Fig.j^a); the triangular area underneath the modal intensity (“IMS triangle”) is 
simply scaled and shifted p9) . 

To understand this observation, we first compare the thresholds (in terms of D 2 ) of both modes 
when evolving the pump profile from different Di values. As we show in Fig.|^b), if we increase 
the D\ value we start with, the non-interacting threshold of mode 1 (in terms of D 2 ) is reduced 
more than that of mode 2, thanks to the stronger intensity of mode 1 in the left cavity. Therefore, 
we expect mode 1 to lase in a wider range of D 2 , which enlarges the IMS triangle. To understand 
the shape invariance of the IMS triangle, we note that it is equivalent to show that the power 
slope of mode I with respect to D 2 does not depend on the fixed Di value along each pump 
trajectory, both before and after the onset of the second mode. 

We resort to SPA-SALT again to explain this property, which also gives the IMS criterion 
when two independent pumps are applied. The constrained two-mode equations 0 now become: 
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where we have used the threshold CF states Q that have the property that one such state is the 
threshold solution, even for a non-uniform index cavity and non-uniform pumping. M here is the 


same as given in 


Eq. 0, 


D, 




is the non-interacting threshold of mode ji with fi = / 2 , and 
1 


W^y = 
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dx (v) M* (x) T7v (x) 


(14) 


is the overlapping factor of mode fi (its corresponding CF state (x) to be exact) with the pump 
in cavity v (v = 1,2), where Tjv(x) = 1 in cavity v and 0 elsewhere in our example above. VF^v 
represents the linear effect we mentioned at the beginning of this section, which indicates that the 
evolving pump profile changes the non-interacting thresholds of the lasing modes. For example. 


D, 


( 1 ), 


( 2 ) 

Q /IF 12 and Dq I'W 21 , whose inverse appear in the matrix P in Eq. (13 1 , are approximations of 


the non-interacting thresholds of mode 1 and 2 when pumping only cavity 2 pT) . The SPA-SALT 
equations dTSll give good approximations of the modal intensity of interacting thresholds as 















Pump strength 

Fig. 5. Invariant shape of the IMS triangle, (a) Same as Fig.^c) but with different values 
of the fixed pump strength D\. Modal intensity of mode 2 is not shown for clarity. Solid 
lines show the full SALT results and dash-dotted lines show the SPA-SALT approximations 
[Eq. f^]. (b) Non-interacting thresholds in terms of £>2 when £>1 is varied. 


Assuming a fixed pump strength D\ and varying D 2 as in Figs.|^c) andj^a), we immediately 
find that the power slope of mode 1 with respect to D 2 , i.e., = fFi 2 is independent 
of the fixed pump strength Di in cavity 1. Likewise, after the onset of mode 2, the power slope 


becomes 

Xii Xn 


(15) 


which is also independent of the fixed Di value along each pump trajectory. This concludes our 
proof of the shape Invariance of the IMS triangle. 

The factor W 22 /Vki 2 in the numerator of Eq. (15 1 makes IMS rely much less on fine mning 
of parameters of the cavity and gain medium. It is basically an indicator that shows how much 
mode 2 (the killer mode here) benefits more from the second pump D 2 than mode 1 (the victim 
mode). For uniform pumping this factor is 1 by definition [see Eqs. Q and ( [T4| )] and we recover 
Eq. Q. In contrast to the uniform pumping case, here even with larger than XUtXii as 

in a typical laser (which gives a positive denominator in Eq. ([fs}), a negative power slope §1 and 
hence IMS can still be Induced by simply focusing the second pump onto mode 2 and avoiding 
mode 1, such that 


1^22 Dphll 
^>2 D^'^xn 


( 16 ) 











As we mentioned previously, X 12 I Xii = 1-64 is larger than = 0.953 in the example 

shown in Figs. and and the criterion ( |10| l for IMS with uniform pumping does not hold. 
However, IF 22 /IF 12 is 2.64 in this case and the criterion ( |16| holds, thanks to the asymmetric 
intensity patterns of mode 1 and 2 in the two cavities [Fig. |^a)]. Equation ( |16[ ) also points 
out the occurrence of IMS by evolving pump profile is more subtle than the intuitive picture 
discussed at the beginning of this section: if the second pump favors mode 2 over mode l(i.e., 
IF 22 /VF 12 > 1) but fails to meet the quantitative requirement imposed by Eq. ( [T^ , IMS does not 
occur either. 

Last but last least, the lasing modes for two identical cavities near threshold are just the 
symmetric and antisymmetric superpositions of the same individual cavity mode in these 
two cavities, and their intensity profiles are both symmetric in the two cavities, resulting in 
W 22 /IF 12 ~ 1 and the absence of IMS, even with independent controls of Di and D 2 . One 
exception is when the system is so lossy that the cavity decay rate becomes larger than the 
coupling between the two cavities, which can result in asymmetric intensity prohles of the two 
modes in these two cavities (see the discussion of the “broken symmetry phase” in p8|). 


5. Conclusion 


We have discussed a mode switching behavior that occurs in lasers with strong non-linear modal 
interactions. In this process the identities of the lasing modes are preserved, including their 
spatial patterns and frequencies. IMS is a robust phenomenon and can in principle occur in any 
type of laser, including standard systems such as a microdisk and coupled one-dimensional 
cavities. In IMS the power slope of a lasing mode turns negative when another lasing mode 
turns on, resulting in the switching off of the first mode in a linear fashion. Qualitatively this 
phenomenon is due to strong cross-saturation effects, and it is explained by the satisfaction of 
two inequalities involving the modal interaction coefficients and the non-interacting thresholds 
of the modes. Specihcally, using the SPA-SALT approximation, we have identihed a criterion 
for IMS with a hxed pump prohle, which requires the cross-interaction coefficient to be stronger 
than the self-interaction coefficient of the killer mode but smaller than that of the victim mode, 
as well as a constraint on their non-interacting thresholds. These conditions can be satisfied, 
for example, when the two modes have similar peak structures and the killer mode is more 
localized. More importantly from the point of view of controlling the mode spectrum of lasers, 
even when these conditions are not satished for a hxed spatial prohle of the pump, IMS can be 
generated by evolving the spatial pattern of the pump as the total pump power increases, no matter 
whether the laser system consists of a single cavity or coupled cavities. This mode switching 
behavior seems to occur in exciton-polariton condensates as well p0| . IMS may have potential 
applications in robust and hexible all-optical switching Eg, stepwise-tunable single-mode laser 
sources, and optical memory cells |32| . An existing experimental setup for its demonstration 
may consist of a microdisk cavity and a spatial light modulator for the optical pump, which has 


been employed successfully to demonstrate pump-controlled emission directionality |331 and 
modal interactions 
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